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I. Introduction: 

In recent years, several authors have investigated the eigenstates of complex potentials Q, es- 
pecially those with PT symmetry and real spectra. In particular, potentials obtained by replacing 
the real coordinate x by a complex variable x + ic in a class of exactly solvable shape invariant 
potentials has been considered by Znojil |2|. For these problems, which are defined on the whole 
real line, an extension to the complex domain was made in order to avoid an inverse-square sin- 
gularity at the origin. The price one pays for taming the singularity is to deal with a complex 
potential. However, it was argued that owing to the PT-symmetric nature of the potential, the 
eigenvalues were still real. As an explicit example, it was shown |^ that a new exactly solvable 
complex harmonic-oscillator like potential with two shifted sets of equally spaced energy levels 
could be generated. The same technique was also applied to explicitly get the eigenstates of 
complex Poschl- Teller I and Poschl- Teller II like potentials Q. 

One of the purposes of this paper is to show that potentials with an inverse square singularity 
at the origin do not necessarily call for moving into the complex domain. In fact, we can obtain 
the spectra of refs. ^ simply by judicious application of the formalism of supersymmetric 
quantum mechanics. Specifically, the spectrum described in ref. Q for the harmonic-oscillator 
like potential is identical to that previously found by us [^, where the discussion focused on 
real potentials with two sets of equally spaced eigenvalues. In this paper, we extend previous 
results and also find new real but singular potentials corresponding to the Poschl- Teller I and 
Poschl- Teller II potentials. Our potentials are shape invariant and consequently their exact 
spectra can be obtained by standard algebraic procedures followed in supersymmetric quantum 
mechanics. We also establish that these singular potentials possess an interesting underlying non- 
linear potential algebra |^, ^ ^, |10| . Explicit representations of the generators are given, which 
provide an alternative algebraic approach to determine the spectrum. 

For completeness, we provide in Sec. 2 a brief review of supersymmetric quantum mechanics 
(SUSYQM)[|l^, 12]. In Sec. 3, we present our framework for generating new shape invariant 



potentials starting from well known solvable problems with an inverse square singularity at the 
origin. We show that if the coefficient of this singularity is restricted within a narrow range, 
one can enlarge the domain of the potential to the negative real axis while maintaining unbroken 
supersymmetry and shape invariance. Working with an explicit example of a harmonic oscillator 
with an inverse square singularity and using the formalism of supersymmetric quantum mechanics. 
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we show that such an extension necessitates the introduction of a 5-function at the origin II which 
eventuahy yields a non-equidistant spectrum for the system. We show that similar extensions 
can be made for Poschl- Teller I and Poschl- Teller II potentials as well and, thus generate new 
shape invariant potentials. We explicitly derive their eigenenergies and eigenfunctions. For these 
potentials, it is important to note that the eigenenergies depend on two parameters, both of 
which get transformed in the shape invariance condition, in contrast to previous work on shape 
invariance. In Sec. 4, we study the potential algebra underlying these systems and generate their 
spectrum by algebraic means. 

2. Supersymmetric Quantum Mechanics: 

In supersymmetric quantum mechanics [|l2[, taking h = 2m = 1, the partner potentials 
V±(x,ao) are related to the superpotential W{x,ao) by 

V±{x,ao) = W\x,ao)±W'{x,ao) , (1) 

where oq is a set of parameters. It is assumed that the superpotential W{x) is continuous and 
differentiable. The corresponding Hamiltonians Hj- have a factorized form 

H = A^A , H+ = AA^ , A = ^ + W{x) , A^ = -^ + W{x). (2) 

ax ax 

We consider the case of unbroken supersymmetry and take Vo ~ exp (— to be nor- 
malizable. This is clearly the nodeless zero energy ground state wave function for H^, since 
AiJo = 0. 

The Hamiltonians and have exactly the same eigenvalues except that has an 
additional zero energy eigenstate. More specifically, the eigenstates of and are related by 

= , Eit\ = 4"^ , oc A V^i-) , A^ V'(+) oc , n = 1, 2, . . . . (3) 

Supersymmetric partner potentials are called shape invariant if they both have the same x- 
dependence upto a change of parameters oi = /(ao) and an additive constant which we denote 
by i?(ao)[^, Often, it is convenient to write this constant in the form of g{ai) — g{ao). The 
shape invariance condition is 

V+{x,ao) = V-{x,ai) + R{ao) = V-{x,ai) + g{ai) - g^ao) . (4) 

The property of shape invariance permits an immediate analytic determination of energy eigen- 
values [|5|, |3|, and eigenfunctions [|^. If the change of parameters oq — > ai does not break 
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supersymmetry, if_(x,ai) also has a zero energy ground state and the corresponding eigenfunc- 
tion is given by -00 0,1) 6xp (^— /^^^ W{y, ai)dyj . Now using eqs. we have 

e(-) =R{ao) , ip[~\x,ao) =A^{x,ao)tp^Q^\x,ao) = ^"^(x, ao) 4"^^;, ai) . (5) 

Thus for an unbroken supersymmetry, the eigenstates of the potential V-{x) are: 

n— 1 n— 1 

=0,4-) = ^ R{ak) = ^ b(afc+i) - giak)] = 9{an) - 5(^0) , (6) 

fc=0 fc=0 

V'i"l\(a;,ai) , (n = 1,2,3, ...) ■ 

These formulas are valid provided the change of parameters oi = /(ao) maintains unbroken super- 
symmetry. In previous work on shape invariant potentials, changes of parameters corresponding 
to translation ai = ao + /? and scaling ai = qao with < q < 1 [^] have been discussed. 
However, a reflection change of parameters oi = —ao, even if it maintained shape invariance, was 
not acceptable since it could not maintain unbroken supersymmetry for the hierarchy of potentials 
built on H_. 

3. New Singular Shape Invariant Potentials: 

The methodology of this paper for obtaining new shape invariant potentials is as follows. One 
begins with a known shape invariant potential, defined for x > 0, which has an inverse square 
singularity A/x^ at the origin. This potential is fully solvable, with eigenfunctions which vanish 
at the origin. One now considers extending the domain to also include the region x < 0. This 
extension is possible only if— l/4<A<3/4. If the strength of the singular term is restricted 
to be in this limited domain, the singularity is called "soft", and the potential is said to be 
"transitional" |T^] . We shall show explicitly how a new change of parameters corresponding to 
the reflection ai = — oq is now admissible, since it maintains both shape invariance and unbroken 
supersymmetry, while still keeping the partner potentials in the soft singularity domain. We can 
then obtain eigenspectra using the shape invariance formalism. As explicit examples, we present 
detailed analyses for the harmonic oscillator, Poschl- Teller I and Poschl- Teller II potentials. 

(a) New shape invariant potential obtained from the harmonic oscillator potential. 

Consider a particle constrained to move in a three dimensional harmonic oscillator potential 

V-{x, 1,00) = ^ uj'^x^ + ^-^^^2^ +(^l-^^uj , {0<x<oo) . (7) 
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-— + T^(x,ao) 
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This potential is generated from the superpotential 

W{x,l,u;) = ^ LOX + - ■ l<0. (8) 

The super symmetric partner potential is 

F+(x,/,a;) = la;V + ^^^+(^^ + ^)a; . (9) 

These two partner potentials are shape invariant since V+{x,l,uj) can be written as V-{x,l — 
l,Lo) + R{l,uj). Here, the remainder R{l,u!) = 2uj is independent of the parameter I. This yields 
an equidistant spectrum En = 2nuj for the harmonic oscillator. However, there is another change 
of parameters that also maintains shape invariance between these two partner potentials, namely 

V+{x,l,Lo) = V-^{x,-l,uj) + R'{l,L0) , 

R'{l,io) = {21 + l)uj. 

However, it is important to point out that since we are at present constrained to be on the 
half-axis x > 0, this second change of parameters, (/,a;) — > is not acceptable for I < 0. 

Neither of the two zero energy solutions -^q"^^ {x,—l,u) oc exp (zt W{x, —I, u) dx) is normalizable 
and hence supersymmetry is spontaneously broken. As shown in sec. 2, the solvability of shape 
invariant systems crucially depends upon supcrpotentials retaining unbroken supersymmetry when 
parameters are transformed, that is, it is essential that V-{x,ai) be a potential with unbroken 
sup er symmetry. 

Let us now consider the same superpotential with an extension of the domain to the entire 
real axis. The asymptotic values of the superpotential are given by the ^ux term at a; — > ±00, 
which is independent of the parameter I. Thus the asymptotic behavior of the ground state wave 
function is dictated by the wx-term and is not affected by flipping of the value of I in the j-term 
of the superpotential. Thus, in contrast with the half-axis case, supersymmetry now remains 
unbroken even with the change of parameters {l,uj) — > (— /, w), and hence this transformation is 
allowed to generate new shape invariant potentials with richer spectra. This leads to 

En = nuj + 2lujPn , P„ = [1 - (-l)"]/2 . (10) 

This new shape invariance yields a new set of eigenenergies superimposed on the old equidistant 
spectrum and are shown in Fig. 1. 
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Figure 1: Energy eigenvalues corresponding to eq. (10) 



We now focus on the region near x = 0. In SUSYQM, it is important that the superpotential 
W{x,aQ) be a continuous and differentiable function. In our example, the above requirement 
is satisfied everywhere except at the point x = 0, where the superpotential of eq. (^) has an 
infinite discontinuity. Such a discontinuity is not acceptable, and needs regularization. Consider 
a regularized, continuous superpotential W{x, ao, e) which reduces to W{x, uq) in the limit e — > 0. 
One such choice is 

W{x,ao,e) = W{x,ao) f{x,e) (11) 



where 



f{x, e) = tanh' 



X 



(12) 



The moderating factor / provides a smooth interpolation through the discontinuity, since it is 
unity everywhere except in a small region of order e around x = 0. In this region, W{x,ao,€) is 
linear with a slope //e^. The potentia]0y_(x, oq, e) corresponding to the superpotential W{x, ao, e) 
is 

V-{x,ao,e) = W\x,ao,e) -W'ix,ao,e) . (13) 



In the limit e ^ 0, y_(x, oq, e) reduces to 



X 



y_(x,ao) = V_{x,ao) -4 Ty(x,ao) —r d{x) , 

X 



(14) 



^At this point one may wonder whether we have lost our cherished shape invariance due to the introduction of 
this moderating factor. In Appendix A, we show that the shape invariance indeed remains intact in the limit e 0, 
and so does the solvability of the model. 
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where we have used lim^^o ^ sech^| = S{x) and linie^o tanh f = yfj- Thus we see that the poten- 
tial V-{x,ao) has an additional singularity at the origin over V-{x,ao) given by ^}{x) = —41 
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Figure 2: The superpotential W of eq. (11) and the corresponding potential V- of eq. (13) for 
the two cases (a) / > and (b) / < 0. 



Note that in the potential shown in Fig. 2(a), the (5-function singularity is instrumental in 
producing a bound state at Eq = 0. 

Naively, in the limit e ^ 0, the potential of eq. (|l^) appears identical to a three dimensional 
oscillator with a frequency lo and angular momentum /. However, there are some more subtle but 
important differences. First, it is defined over the entire real axis (— oo < x < oo) and not just 
the half line. For a proper communication between the two halves, we must have a "softness" of 



the inverse square term. Normalizability of the wave function requires that the coefficient A of the 

inverse square term be in the transition region — ^ |jl^. More specifrcahy, for > 0), one has 

< /(/ + 1) < I and for (/ < 0) one has < l{l + 1) < 0. The important special case of the one 
dimensional harmonic oscillator has / = 0: it corresponds to l{l + 1) = and no singularity. 
For transition potentials, both solutions of the Schrodinger equation are square integrable at the 
origin. Therefore, both are acceptable square integrable solutions of the Schrodinger equation, 
and must be retained to form a complete set. Eigenstates for the potential V-{x,aQ) can be 
obtained from eq. (^. The lowest four are 



^0 = 0; V'o oc x"' , 

El = {21 + V^i x^+^ e"^"^^'" , 

E2 = 2uj ; ^2 oc f 2/ - 1 + ux'^) x-'e"^"^"^" 



E3 = 2uj + (21 + l)uj ; V'a (-2/ - 3 + 



(15) 



General expressions for these eigenfunctions and corresponding eigenenergies are 

„2' 



E2n = 2nu! , 



V'2n OC X^'e"^'^^ Ln ^ 



LOX 



E2n+i = 2nu + {21 + 1)W , V2n+1 OC x^+'c" 3"^^' L^"^^ 



where L„ are the standard Laguerre polynomials. 



uox 



(16) 



(b) New shape invariant potential obtained from the Poschl- Teller I potential. 

As a second example, we consider the Poschl- Teller I superpotential 

W{x,A,B) = Aia.nx-Bcoix ; < x <it/2 . (17) 
The super symmetric partner potentials are then given by 

V-{x,A,B) = A{A- 1)860"^ x + B{B-l) cosec^x- {A + Bf , (18) 

and 

V+{x,A,B) = A{A + l)sec'^x + B{B + l)cosec'^x-{A + Bf . (19) 
Here, A and B are both positive in order for V-{x, A, B) to have a zero energy ground state. 



Again, one can readily check that there are two possible relations between parameters such that 
above two potentials exhibit shape invariance. One of them is the conventional {A ^ A + 1, B ^ 
B + 1). The second possibility \s {A ^ A + \, B ^ ~B). As explained in the previous section, 
this second relationship breaks supersymmetry on (0, 7r/2) domain and it is allowed only if the 
domain of x is extended to range — 7r/2 < x < 7r/2. The first transformation among parameters 
{A ^ A + \, B ^ B + 1) has been studied extensively in the literature. It is the second 
transformation that yields new results and will be considered here. Thus, the relationship among 
parameters that we consider is, (^fc+i = Ak + 1, -Bfc+i = —Bk). This potential also requires a 
careful analysis in the vicinity of x = 0, where two half-axes are being sewed together. Again, 
the need of continuity and differentiability of the superpotential requires its regularization, as was 
done in eq. for the harmonic oscillator. A similar analysis then leads to a new singular shape 
invariant potential 

V- ix, A, B) = \A{A - 1) sec^ x + B{B-l) cosec^x -{A + Bf] + 4 5 cot x Six) . (20) 

1^1 

This potential obeys the shape invariance condition: 

V+{x,A,B)=V.ix,A + l,-B) + {A + l-Bf-{A + Bf , (21) 
and its eigenvalues and eigenfunctions are given by 
Eo = 0, 

tpQ (X COS X sm X , 

Ei = {A + 1- Bf - {A + Bf , 

tpi oc cos"^ X sin~'^~^ x [{2B — 1) cos^ x + 1] , 

E2 = {A + 2 + Bf - {A + Bf , 

ip2 oc cos^ X sin-^ x [{AB + 2) cos^ x - (65 + 3) cos^ x + 1] , 
E3 = {A + 3- Bf - {A + Bf , 

tp3 OC cos^ X sin"-^"^ x U-SB"^ + 16B - 6) cos*^ x + (125^ - 32B + 13) cos^ x + {20B - 8) cos^ x + 1 
Thus, the general formula for eigenvalues is 

Er, = {A + n+i-irBf-iA + Bf . (22) 
The eigenspectrum is shown in Fig. 3 and Fig. 4. 







1 






1 








-\- 




4— 


- 




1- 


1 




■ 






















\ 1 

















Figure 3: The potential of eq. (pO|) for A = 1.5 and = —1/3 and its energy spectrum. 

Note that, to avoid level crossing, we must have En > En-i- This leads to the constraint 
— ^ < i? < ^. Interestingly, it is the same constraint that one needs for the normalizability of 
the wavefunction at the origin and hence to the possibility of communication between regions 
(— |7r,0) and (0, ^vr) of the domain. 



(c) New shape invariant potential obtained from the Poschl- Teller II potential. 

The last example that we consider is that of the Poschl- Teller II potential described by 

W{x,A,B) = Atanh.x-Bcot'hx ; < x < oo . 



(23) 



Here, A and B both need to be positive and satisfy the condition A > B for the potential 
V-{x, A, B) to have a zero energy ground state and to ensure unbroken supersymmetry. The 
supersymmetric partner potentials are then given by 



V+{x, A, B) = -A{A - 1) sech^x + B{B + 1) cosec^x + {A - Bf 



(24) 



and 



V_{x, A, B) = - A{A + I) sec^ x + B{B- I) cosec^x + {A -Bf . (25) 

Here too we have two possible relations between parameters for these potentials to be shape 
invariant. They are (A — > ^4 — 1, i? — > S + 1), and (^ ^ ^ — 1, B ^ ~B). As explained before 
in last two examples, the second transformation requires an extension of the range to (— oo, oo). 
The new singular potential generated for this case is given by 



V-{x,A,B) 



-A(A - 1) sech^x + BiB + 1) cosech^x + (A - Bf +4B coth x — 6{x) . (26) 

J X 
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Figure 4: The superpotential W potential V- from Poschl-Teller for A = 1.5 and B = —1/3. 
The shape invariance condition obeyed by this potential is given by 

V+{x,A,B) = V-{x,A + l,-B) + {A-Bf-{A-l + Bf , (27) 
and the eigenvalues and eigenfunctions are given by 
Eo = 0, 

tpo oc cosh""^ X sinh^ x , 

Ei = {A-Bf -{A-1 + Bf , 

ipi oc cosh~^ X sinh~^~^ x [{2B - 1) cosh^ x + 1] , 

E2 = {A- Bf - {A-2- Bf , 

'02 OC cosh"^ x sinh^ x [(4i? + 2) cosh^ x — {QB + 3) cosh^ x + 1] , 
E3 = {A-Bf-{A-3 + Bf , 



V's oc cosh^ X sinh"^-^ x {8B^ - 16B + 6) cosh^ x - (125^ - 325 + 13) cosh'' x 

-(20B-8)cosh2x- ij . (28) 

Thus, the general formula for eigenvalues is 

En = {A-Bf -{A-n-{-l)''Bf . (29) 

Again, to steer clear of the level crossing problem, we must have En > En~i- This leads to the 
constraint — ^ < i3 < ^; which, as stated earlier, is the same constraint that one needs for the 
normalizability of the wavefunction at the origin and for an effective communication between two 
halves of the x-axis. 

4. Potential Algebra: 

So far, we have discussed three types of new solvable singular potentials. We will now derive 
the potential algebra underlying them. We will show that the algebra based on the generators 
{ J+, J_, J3} is non- linear [^, @, H], 10 1 . Potential algebras provide an alternative way of getting 



the eigenvalues by algebraic means. 
Consider the following ansatz: 

J^ = c-^ AHx,a{N),P{N)) , J.=A{x,a{N),f5{N))c, J^ = N = c^c, (30) 

where c, and are three operators satisfying [c, c^] = 1, and cc^^=c^^c = l. An example 
of such operators is given by c = e*"^, = e"*"^ and = iS^e"*'^, where 4> is some arbitrary 
real variable. The operators A and A"^ of eq. ( 30) are obtained from eq. @i via the substitution 
oq = {A, 5} {a(A^), /3(A^)}, where a and [3 are real, arbitrary functions to be determined later. 
We can readily check that 

[J3,J±] = ±J±, [J+,J-] = -R{,h)=g{a{N),P{N))-g{a{N -l),p{N -I)) . (31) 

The last commutation relation is a consequence of the algebraic shape invariance condition 

H+{x,a{N),P{N))-H^{x,a{N -1),(3{N -I)) = g{a{N -1),(5{N -l))-g{a{N),(5{N)) , (32) 

which is the operatorial "twin" of the classical shape invariance condition eq. (^ obtained via 
the mappings {Aq.Bq} {a{N),(3{N)} and respectively {Ai,Bi} {a{N - 1),(5{N - 1)}. 

The functions a{N) and (3{N) are determined by requiring that the change a{N) a{N — 
1) and [3{N) [3{N — 1) correspond to the change of parameters oq — > ai. For example, 
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a{N) = A — N corresponds to a translational change of parameters Aq ^ Ai = Aq + 1, because 
a{N — 1) = a{N) + 1. Similarly, f3{N) = (—1)^ corresponds to the reflection Bq —>■ Bi = —Bq , 
since (3{N - 1) = -/3{N). 

For any shape invariant potential, we know the function g(a,l3), which explicitly gives the 
potential algebra (^Tj). From its representations, we can obtain the energy spectrum for the given 
problem. 

To find a representation of the potential algebra, let us consider a set of eigenvectors common 
to both = J-i- J_ and J3 = N denoted by {|n), n = 0, 1, . . .}. The action of J+, J_ and J3 on 
this basis is given by 

J+|n) = a(n + l)|n + 1) , J-|n) = a(n)|n — 1) , Jsjn) = n|n) . (33) 

Here we have chosen, without any loss of generality, the coefficients a{n) to be real. Note that 
since J-|0) = 0, we have the initial condition a(0) = 0. There is a connection between the 
coefficients a(n) and the eigenspectrum of the Hamiltonian. Observe that 

H_{x,a{N - l),/3(iV - 1)) \n) = J+J_|n) = a^{n)\n) . (34) 

Therefore, in order to find the spectrum of the Hamiltonian we have to determine the coefficients 
a^(n). This can be done by projecting the last equation from (^Tj) on \n) and solving the resulting 
equation recursively. Thus, we obtain a-^(n)— a^(n+l) = g{n)—g{n—l) having the solution a?{n) = 
g{—l) — g{n — 1). Here we have denoted g{n) = g{a{n), (3(n)). But a?{n) corresponds to the 
eigenvalues of the Hamiltonian H-{x, a{N — 1), (3{N — 1)), or "classically" speaking to the shifted 
set parameters oi. Therefore the eigenenergies of the initial Hamiltonian H^{x,a{N), (i{N)) 
(corresponding to the set of parameters gq) are 

En = g{a{0),m) - 9{ain),f3{n)) . (35) 
We make contact with eq. (^) by observing that {a{n — k),(3{n — A:)} = a^. 

(a) The harmonic oscillator. 

To show how our procedure works, it is instructive to build explicitly the potential algebra 
of the harmonic oscillator. The superpartner potentials and V+ are given in eqs. (0) and 
respectively (^). Under the change of parameters {/,a;} — > {/ — we have the following shape 
invariance condition 

H+{x,l,u;) = H^{x,l- l,u;) + {-2lo {I - 1)) - {-2ujI) . 
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To build the potential algebra, first we find the functions a and /3 associated with the above 
change of parameters. We have immediately a(A^) = / + , /?(A^) = Next, we can build the 
concrete realization of the potential algebra using the ansatz ( pO[ ) and the superpotential from 
(|8|). The resulting generators 

1 / d 1 l^N\ ^ f d 1 l + N\ , + 

J, = c-(-- + -.. + ^j , ^-=(^ + r^ + -^j^' ^3 = A^-ctc, (36) 

satisfy the "canonical" commutation relations ([3l|), where the function g is given by g{N) = 



g{a{N), ()(N)) = —2uj [l + N). Finally, using the formula (35) we get the spectrum £"„ = g{Q) — 
g{n) = 2a; n, which is exactly what we have expected. 

Next, let us consider the new singular shape invariant potential corresponding to the change 
of parameters {/,c<;} {— Z,c<;}. In this case a{N) = —{—1)^1 and (3{N) = u). From eq. (^) we 
get 

^ ( d 1 , ( d I (-l)^A , + , 

J, = c- + '^-=(^ + r^-^j^' ^3 = iV^ctc.(37) 

The commutation relations (^) together with the algebraic shape invariance condition (^2|) yield 
in this case [J-|_, J_] = —uj (— 2(— 1)^/ + 1) from where we get ^(A^) = w ((—1)^/ — N). Therefore, 
the resulting eigenspectrum (pBf) is = ujn + — (—1)"). 



(b) The Poschl- Teller I potential. 

We build the algebraic model for the new shape invariant Poschl- Teller I like potential by 
taking into account that corresponding to the change of parameters {A,B} — > {A + 1, — -B} we 
have a{N) = ^ - and (3{N) = {-l)^B. Then, using the superpotential (0) one gets the 
following expressions for the generators of the associate potential algebra 

J+ = c-^ + (yl - iV) tan x- (-1)^5 cot , 

j_=(-^j^{^A-N)tSinx-{-l)^Bcoix\ c, J3 = N = c^c. (38) 



Using as before the algebraic shape invariance condition (32) we obtain in this case [J+, J„] = 
-{A + N + l + {-l)^^+^^Bf + {A+N + {-l)^Bf. Therefore we get g{N) = -{A + N + Bf 
and the corresponding eigenspectrum £'„ = g{0) — g{n) = —{A + B)"^ + {A — n + {—1)^B). 

(c) The Poschl- Teller II potential. 
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For the new the Poschl- Teller II potential like case, to the change of parameters {A, B} — s- 
{A — 1, —B} we have a{N) = —{A — N) and (3{N) = {—1)^B and the corresponding algebra is 
therefore generated by 

J+ = c"^ (^--^ + {-A + N)tanhx - i-l)'^Bcothx^ , 

J_ = f--^ + {-A + N)tauhx - i-l)^Bcothx] c , J3 = N = c^c. (39) 



In the above representation the explicit form of the superpotential ( |2^ ) was taken into account. 
The commutation relations (^) together with the algebraic shape invariance condition (^2[) yield 
in this case ^(A^) = {—A + N — (— Using (pS]), one obtains as expected, the eigenspectrum 
for this potential = 5(0) - g{n) = {A + Bf - {-A + n- {-l^Bf. 

5. Conclusions and Comments: 

We have generated several new shape invariant potentials on the whole line starting from well 
known potentials on the half line. To ensure continuity and differentiability of the superpotential, 
our procedure requires a regularisation at the origin. This extension not only maintains shape 
invariance, it also allows the possibility of a new transformation among parameters {B — > —B) 
that was not allowed on the half-axis. This transformation results in new superpotentials, albeit 
singular, that are defined over the entire real axis and have richer spectra than those defined over 
half-axis. It is shown further that the eigenspectra of these new real singular shape invariant 
potentials may also be derived from a nonlinear potential algebra. 

Since we have obtained and discussed the exact eigenvalues and eigenfunctions of three new 

singular potentials using the machinery of supersymmetric quantum mechanics, it is of interest 

to ask what one gets in the WKB approximation. Let us recall that Comtet et al. have shown 

the exactness of the SWKB quantization condition[18, 13] 

rx2 

/ yjEn -W^ dx = mrh . 

J XI 

for all known shape invariant problems with unbroken SUSY where parameters are related by 



ai = ao + 6 [21 



For broken SUSY, Inomata and Junker gave the quantization condition 

n VEn -W^ dx=[n + l/2]7Tn . 

J XI 

For both cases, the turning points xi,X2 are solutions of W'^{x) = En- 
lb 



Our new singular potentials allow a change of parameters that, if considered in half- axes 
only, leads the system to alternate through unbroken and broken phases of supersymmetry as 
o-k fflfc+i- It is interesting to note that the spectrum of these singular potentials can be derived, 
using somewhat more complex but exact quantization condition which alternates between the 
broken and unbroken SUSY cases: 

/ VEn -W^ dx = [n + 1/2P„] , 

where P„ is equal to [1 - (-l)"]/2. 

Partial financial support from the U.S. Department of Energy is gratefully acknowledged. 
R.D. would like to thank the Department of Atomic Energy, Government of India for a research 
grant, and the Physics Department of the University of Illinois at Chicago for warm hospitality. 



Appendix A: In this appendix, we show that shape invariance of our new potentials is main- 
tained during the process of extending the domain to the whole real axis and introducing the 
moderating factor f{x,e). Let us recall that our old superpotential W{x,A,B) is of the form 
A -|- B [^(x)]^^, where the function is x,tana; or tanhx for harmonic oscillator, 

Poschl- Teller I, and Poschl- Teller II respectively^ 

Note that in all cases, ^ — >l/x at the origin. W{x, A, B) is replaced by a regularized, 

continuous superpotential W{x^A^B,e) given by 

W{x, A, B, e) = W{x, A, B) f{x, e) , (41) 

where f{x,€) is unity everywhere except in a small region of order e around x = 0. One such 
function f{x,e) is given by tanh^ i^/^)- Iii limit e — > 0, we assume that the f{x,e) 1 and 
'^•^^^'^^ — > 2-^6{x). The potentials V^{x, A, B, e) corresponding to the superpotential W{x, A, B, e) 
are then given by 

V^{x, A, B, e) = W\x, A, B)f\x, e) T ( "^^^^^'Z' /(^, ^) + W{x, A, B)) . (42) 



^The change of parameters associated with shape invariance in these potentials are of the form 



and B — y B = -B . ■ (40) 
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Now 



V+(x,A,B,e)-V-(x,A,B,e) = V+(x,A,B)-V-{x,A,B) + ^^^^ (w(x,A,B) + W(x,A,B)) 

ax ^ ) 

= R(A,B) + ^i^ {B + B)Mxr' 

= R{A,B) , (43) 

where we have used the hmits of / and /' and B = —B. This estabhshes the shape invariance of 
the regularized superpotential. 
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